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Generic model predictions for the early neutrino signal of a core-collapse supernova (SN) imply 
that IceCube can reconstruct the bounce to within about ±3.5 ms at 95% CL (assumed SN distance 
10 kpc), relevant for coincidence with gravitational- wave detectors. The timing uncertainty scales 
approximately with distance-squared. The offset between true and reconstructed bounce time of 
up to several ms depends on the neutrino flavor oscillation scenario. Our work extends the recent 
study of Pagliaroli et al. [Phys. Rev. Lett. 103, 031102 (2009)] and demonstrates IceCube's superb 
timing capabilities for neutrinos from the next nearby SN. 
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I. INTRODUCTION 

The high-statistics neutrino observation from the next 
nearby supernova (SN) will provide a bonanza of infor- 
mation about the astrophysics of core-collapse phenom- 
ena and neutrino properties. In a recent Physical Review 
Letter a strong case was made for the importance of co- 
incidence measurements between gravitational wave and 
neutrino signals from SN core bounce [l[ . The largest ex- 
isting SN neutrino detectors arc Super-Kamiokande and 
IceCube, reaching to a distance of about 100 kpc. The 
expected distribution of galactic SNe drops off quickly 
beyond about 20 kpc At this "pessimistic" distance, 
Super-Kamiokande can time the bounce to within a few 
tens of milliseconds, an interval comparable to the ex- 
pected duration of the gravitational wave burst [l| . 

We here extend this study to IceCube, a high-statistics 
SN neutrino detector that would have seen the SN 1987A 
neutrino signal with 5 a significance. For galactic SNe, 
the large rate of uncorrelated Cherenkov photons pro- 
vides excellent time-structure information. 

Following Ref. [H we note that neutrino masses are 
small enough to neglect time-of-flight effects. Recalling 
that the diameter of the Earth is 42 ms, millisecond- 
coincidence measurements between detectors at differ- 
ent geographic locations requires determining the SN 
direction either by astronomical observations or by the 
electron-scattering signal in Super-Kamiokande. 



in 1.6384 ms bins, implying a total event rate of 

R^^ - 186 bin^^ L52 D^^ {El^)/{Eii) , (1) 

where L52 = Lp_^/1G^^ erg s^^, fio = D/IQ kpc and 
EiT, = -EpJ15 MeV. 

The peak luminosity reaches L52 = 2-5 and at that 
time {EjyJ) w 15 MeV. For a thermal or slightly pinched 
spectrum, (i?? )/{Ep^)^ « 2. Altogether, we expect 
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1.5 X 10-^ bin" 
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as a typical peak event rate for a SN at 10 kpc. 

This signal is to be compared with a background of 
280 s-i in each OM Q, corresponding for 4800 OMs to 



Rn = 2.20 X lO'^ bin" 



(3) 



with an rms fluctuation of 47 bin" . Therefore, the SN- 
induced "correlated noise" in the entire detector is highly 
significant 0, @, @] . From any one neutrino interacting in 
the ice, at most one Cherenkov photon is picked up: the 
signal counts are entirely uncorrelated. 

The prompt Ve burst immediately after bounce pro- 
duces a peak rate of about 100 bin~^, lasting for several 
ms. Its impact on the early count rate depends sensi- 
tively on the flavor oscillation scenario (see below). 



III. EARLY NEUTRINO EMISSION 



II. SUPERNOVA NEUTRINOS IN ICECUBE 

When SN neutrinos stream through water or ice, 
Cherenkov light is generated, primarily by the secondary 
positrons from inverse beta decay i^e+P — > n-\-e'^ . While 
fine-grained detectors reconstruct individual neutrinos on 
an event-by-event basis, IceCube only picks up the aver- 
age Cherenkov glow of the ice. To estimate the detection 
rate we follow Ref. Q, augmented with the latest Ice- 
Cube efficiencies [1]. The complete detector will have 
4800 optical modules (OMs) and the data are read out 



SN models suggest that neutrino emission for the first 
20 ms after core bounce depends little on model assump- 
tions or input physics Q, although beyond this early 
phase the accretion rate and therefore neutrino emission 
depends strongly, for example, on the progenitor mass 
profile. The prompt v^. burst reaches its peak at 5-7 ms 
post bounce, which also marks the onset of De emission 
that is initially suppressed by the large chemical po- 
tential. Up to about 20 ms post bounce rises roughly 
linearly. Wc thus represent the early IceCube signal as [l| 

n pmax ^ / for t <. , .s 

X 1 1 _ e-(*-*')/-.- fort>t, 
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with U = 6 ms, t, = 50 ms and Rf^'''' = 1.5 x 10^ bin"^ 
These parameters also provide an excellent fit to the first 
100 ms of a numerical model from the Garching group Q 
that is available to us. 

We may compare these assumptions with the early- 
phase models of Rcf. Q- L^^ rises nearly linearly to 
L52 = 1.5-2 within 10 ms. The evolution of (i?p^)RMS = 
{{E^J / {Ep^})^/"^ is also shown, a common quantity in 
SN physics that characterizes, for example, the efficiency 
of energy deposition; the IceCube rate is proportional 
to (£'pe)RMS- 10 ™s after onset, (i?p,)RMS reaches 
15 MeV, implying {Ef^)/{Ei5) = 1. We thus estimate 
10 ms after onset a rate of 280-370 bin~^, to be compared 
with 270 bin^^ from Eq. (|4]). Therefore, our assumed sig- 
nal rise is on the conservative side. 

Of course, the early models do not fix t^- and i?™^^ 
separately; the crucial parameters are and R^'^^/t^. 
The maximum rate that is reached long after bounce is 
not relevant for determining the onset of the signal. 

If flavor oscillations swap the flux with Dx (some 
combination of and Vt), the rise begins earlier be- 
cause the large chemical potential during the pronipt 
Uf. burst does not suppress the early emission oi Q. 
Moreover, the rise time is faster, (£')rms larger, and the 
maximum luminosity smaller. We use Eq. (|4]) also for 
with t, = 0, T, = 25 ms, and Rf^"" = 1.0 x 10^ bin"\ 

Flavor oscillations are unavoidable and have been stud- 
ied, for early neutrino emission, in Ref. 0. Assuming 
the normal mass hierarchy, sin^ 613 ^ 10""^, no collec- 
tive oscillations, "'^ and a direct observation without Earth 
effects. Table I of Ref. Q reveals that the burst would 
be completely swapped and thus nearly invisible because 
the 1^1:6" elastic scattering cross section is much smaller 
than that of Ve- The survival probability of D,, would be 
cos^ 012 w 2/3 with 812 the "solar" mixing angle. There- 
fore, the effective detection rate would be | Riy^ + ^ . 
We use this case as our main example. 



IV. RECONSTRUCTING THE SIGNAL ONSET 

A typical Monte Carlo realization of the IccCubc signal 
for our example is shown in Fig. [1] One can determine 
the signal onset to within a few ms by naked eye. For a 
SN closer than our standard distance of 10 kpc, one can 
follow details of the neutrino light curve without any fit. 

One can not separate the and components for 
the example of Fig. [T] Therefore, we reconstruct a fit 
with a single component of the form Eq. ([4]), assuming 
the zero-signal background is well known and not fitted 



^ In the normal hierarchy, collective oscillation effects are usually 
absent. It has not been studied, however, if the early neutrino 
signal can produce multiple splits that can arise also in the nor- 
mal hierarchy Moreover, for a low- mass progenitor collective 
phenomena can be important if the MSW resonances occur close 
to the neutrino sphere llOl [TUl . 




Time post bounce (ms) 

FIG. 1: Typical Monte Carlo realization (red histogram) and 
reconstructed fit (blue line) for the benclimark case discussed 
in the text for a SN at 10 kpc. 



here. Using a time interval until 100 ms post bounce, 
we reconstruct = 3.2 ± 1.0 ms (Icr). If we use only 
data until 33 ms post bounce we find <o = 3.0 ±1.7 ms. 
Indeed, if one fits Eq. @ on an interval that ends long 
before the plateau is reached, we effectively fit a second 
order polynomial with a positive slope and negative sec- 
ond derivative at fr, whereas the plateau itself is poorly 
fitted and its assumed value plays little role. Depending 
on the distance of the SN one will fit more or fewer details 
of the overall neutrino light curve and there may be more 
efficient estimators for t-^. Our example only provides a 
rough impression of what IccCube can do. 

The reconstruction uncertainty of to scales approxi- 
mately with neutrino flux, i.e., with SN distance squared. 
The number of excess events above background marking 
the onset of the signal has to be compared with the back- 
ground fluctuations. Therefore, a signiflcant number of 
excess events above background requires a longer integra- 
tion period if the flux is smaller, explaining this scaling 
behavior. 

The interpretation of <o relative to the true bounce 
time depends on the flavor oscillation scenario realized in 
nature. This is influenced by many factors: The value of 
Oi3, the mass ordering, the role of collective oscillation 
effects, and the distance traveled in the Earth. Com- 
bining the signal from different detectors, using future 
laboratory information on neutrino parameters, and per- 
haps the very coincidence with a gravitational-wave sig- 
nal may allow one to disentangle some of these features. 
However, as a first rough estimate it is sufficient to say 
that the reconstructed to tends to be systematically de- 
layed relative to the bounce time by no more than a few 
ms. The statistical uncertainty of the t^ reconstruction 
docs not depend strongly on the oscillation scenario. 
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V. DISCUSSION 

The authors of Ref. [l| used a SN emission model based 
on a two-component fit of the sparse SN 1987A data and 
described the overall neutrino signal in terms of several 
parameters which they say, after Eq. (5), are at odds 
with theoretical expectations. However, their parame- 
ter Ta = 2.4 MeV, describing the temperature of the 
neutrino-emitting gas during the accretion phase, gives 
(-Bp J = 5.2 = 12.5 MeV (seeparagraph after Eq. 15 
in the published version of Ref. [12]) and thus is virtually 
identical to the corresponding {EpJ^iMS from Ref. (3] at 
20 ms post bounce. In other words, while the SN 1987A 
implied z/g energies of Ref. [l[ are lower than theoreti- 
cal expectations for the overall accretion phase, in the 
absence of flavor oscillations they agree nicely with the 
models of Ref. 0] for the crucial first 20 ms. However, 
the chosen rise time Tr = 100 ms is very long compared 
with the early-time models of Ref. . 

In Ref. [3l the early neutrino signal after bounce was 
systematically studied for different input assumptions 
(progenitor mass, equation of state, neutrino opacities), 
leading to very similar results. On the other hand, 
one finds significantly different numerical examples in 
the literature. In Ref. [l^ the Pe luminosity rises to 
2 X 10^^ erg s~^ after as much as 50 ms. In Ref. (l3 |. 
a peak value of only 1 x 10^^ erg s~^ is reached and the 



rise within 10-20 ms after bounce is small. In Ref. [15 1 
the signal begins rising as late as 12 ms after the maxi- 
mum of the Ve burst. There are many differences between 
these and other models in terms of numerical approach 
and input physics. It would be extremely useful if an- 
other group would investigate the early neutrino signal 
in the spirit of Ref. (3] for a range of physical assump- 
tions and with attention to numerical details that may 
influence the early-time behavior. 

In view of the large range of possible distances to the 
next nearby SN and concomitant flux differences, these 
uncertainties do not change our overall conclusions. In 
our fiducial example IceCube can reconstruct the signal 
onset within ±6-7 ms at la CL for a SN at 20 kpc, com- 
parable to what Ref. [l[ found for Super-Kamiokande. 
Ideally, of course, one would combine the measurements 
from several detectors. 

A gravitational wave measurement of the core bounce 
in coincidence with neutrino onset would be of obvious 



astrophysical importance. In addition one could test the 
weak equivalence principle. Both neutrinos and gravita- 
tional waves should suffer the same Shapiro time delay in 
the gravitational potential of the galaxy. For SN 1987A in 
the Large Magellanic Cloud, this delay was a few months. 
The coincidence of the neutrino burst with the rise of the 
light curve within a few hours proved an equal Shapiro 
delay for photons and neutrinos to within about 10"'^ 
16l . 17 1 . A millisecond-scale coincidence between neutri- 
nos and gravitational waves would extend and refine this 
test, in detail depending on the location of the SN. 

If the SN can not be located by astronomical means be- 
cause of obscuration, the electron-recoil signal in Super- 
Kamiokande or a future megatonne water Cherenkov de- 
tector is the method of choice [3, [l^ , whereas arrival- 
time triangulation was dismissed. With IceCube al- 
most complete, the situation has changed and triangu- 
lation can play a useful role after all, at least for a 
not-too-distant SN. Moreover, in Europe a big detec- 
tor of the LAGUNA-class [20] may become available 
(representing one of three possible large-scale detector 
types) and a megatonne water-Cherenkov detector may 
be built in North America. Arrival-time triangulation 
between large-scale neutrino detectors on different conti- 
nents could become viable. 

It would be a worthwhile exercise to study the possi- 
bilities of bounce timing, based on common astrophysical 
assumptions, for different detector types and what can be 
learnt from a combined analysis. To this end it would be 
worthwhile if groups other than the Garching one would 
systematically study the numerical early-time neutrino 
signal to judge if indeed the dependence on physical as- 
sumptions is as small as found in Ref. [3]- 
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